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SUMMARY 


An investigation vas made to determine the effect of longitudinal- 
stiffeners on the buckling load of long flat plates under shear. 

Tests were made of long flat plates reinforced by one and by two 
longitudinal stiffeners. A theoretical study of the buckling load of 
such plates, made by the energy method, is presented in the appendix. 
The results of the tests and the results of the theory are compared 
and are found to be in fair agreement. 


INTRODUCTION 


Buckling of the stressed skin of a wing under applied shear 
loads results in a reduced torsional stiffness and a reduced aero- 
dynamic fairness of the wing. Because there is danger of flutter or 
aileron reversal occurring if the torsional stiffness of the wing is 
not maintained up to high loads and because high-speed pull-outs may 
be difficult or impossible if reduced aerodynamic fairness causes 
premature separation of the flow over the wing, it is desirable to 
determine the shear stress at which the reinforced skin of the wing 
buckles. The problem is of particular importance in the case of 
high-speed airplanes which are normally subject to flutter and 
control problems. With a view toward eliminating some of the 
problems in high-speed flight, therefore, a solution to the problem 
of the shear buckling of a type of panel likely to be used in the wings 
of fast airplanes has been sought . 

The thin wings needed for high-speed airplanes have thick skins 
and several shear webs. The wing panels are narrow and, therefore, 
are reinforced by relatively few stiffeners. Accordingly, tests were 
made to determine the shear buckling load of long plates reinforced 
by one and by two longitudinal stiffeners. In addition, a theoretical 
solution of the problem for any number of stiffeners was made. The 
results of the tests are presented herein and are compared with the 
results of the theory. 

The symbols are defined in the appendix. 
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TEST SPECIMENS 


The specimens tested were flat plates with a length— width ratio 
of 8 reinforced by longitudinal stiffeners. The general construction of 
the test specimens is shown in figure 1, and the specific dimensions of 
the individual specimens are listed in table I. 

Two groups of specimens were tested. Each specimen of the first 
group had one stiffener riveted along the longitudinal center line of the 
panel. Each specimen of the second group had two longitudinal stiffeners 
of equal size riveted to the panel in order to divide the panel into three 
hays of equal width. The dimensions of the specimens of the first group 
were nominally 6 inches wide and 48 inches long. For the dimensions of 



the attached legs of the stiffeners would not cover a large part of the 
width of the hays, and the length was increased to 63 inches in order that, 
at the same time, the length-width ratio should remain at a value of 8. 

The wehs of all specimens were of nominally O.O32— inch— thick 
24S— T aluminum— alloy sheet. The stiffeners were of 24&-T aluminum— alloy 
sheet hent to the shape of angles. Two angles were used for each stiff- 
ener, one on each side of the weh, to provide symmetry about the plane of 
the weh. The thickness and leg dimensions of the angles were varied to 
produce the bending stiffness desired. 

The short edges (ends) of the specimens were reinforced with angles. 
These angles were of uniform size for all specimens with one stiffener and 
were proportionally larger and of uniform size for all specimens with two 
stiffeners. These end angles were so designed that there was a margin of 
safety against failure of the angles before buckling of the web occurred. 


TEST APPARATUS AND TESTING PROCEDURE 


The specimens were tested in a Jig as shown in figure 2. One part of 
this Jig distributed the applied load along one edge of the web, and the 
other part picked up the reaction from the opposite edge of the web and 
transferred this reaction to a heavy supporting structure. Both parts 
of the jig were essentially the same, and each part consisted of two 
heavy steel bars bolted to each side of a steel plate which protruded 
from between the bars and to which the specimen was riveted. The large 
cross-sectional area of the bars insured that the distribution of load 
was essentially uniform over the full length of the web. 

A portable hydraulic jack which indicated loads with standard 
testing-machine accuracy of one— half of 1 percent, was used to apply the 
load to the specimens. 

Two dial gages graduated to l/lOOOO inch were used to measure the 
shear displacement of the loaded edge of the sheet relative to the fixed 
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edge. These gages were mounted on each side of the wet at the midlength 
of the specimen. 

The test procedure was as follows: The specimen was preloaded in 

several increments to about 25 percent of the estimated buckling load, and 
the dial gages were read after each increment of load had been applied. 

If the dial— gage readings indicated equal movement on both sides of the 
web up to the full preload, it was assumed that the jack was properly 
positioned under the specimen and the load was released. The load was 
reapplied and the dial— gage readings were taken at a number of loads until 
buckles were clearly visible in the web. Leadings were then taken less 
frequently until the specimen could sustain no further increase in load. 


ANALYSIS AND DISCUSSION 


Theoretical critical stresses .— The theoretical study of the shear 
buckling load of a long flat plate , reinforced by longitudinal stiffeners, 
with shear load acting at the longitudinal edges is presented in the 
appendix. Two different restraining conditions at the longitudinal edges 
are investigated — simply supported edges and clamped edges. The results 
of the theoretical study are summarized in figure 3 . For each edge 
restraint condition, three separate curves are shown. The curves corres- 
pond to the plate reinforced by one stiffener, by two identical equally 
spaced stiffeners, and by a large number of identical equally spaced 
stiffeners. The three curves are essentially one with a maximum deviation 
of approximately 2 or 3 percent. The appendix points out that the curve 
for a large number of identical equally spaced stiffeners can be used to 
represent, within 2 or 3 percent, the solution of the plate reinforced by 
any number of identical equally spaced stiffeners. 

Each curve in figure 3 has a * 1 upper limit corresponding to the crite- 
rion that the plate buckles in such a way that the stiffeners can be 
replaced by simple supports. In the case of a simply supported plate with 
one stiffener, for instance, the shear buckling coefficient k cannot be 
increased beyond approximately 21.4 by increasing the stiffness of the 
stiffener. 

Expe rimenta l buckling data.— In figure 4 are shown typical results of 
a test in the form of a curve of shear deformation, as measured by the 
dial gages, plotted against load on the specimen. The first part of this 
curve is linear (i.e., deformation is proportional to load) and corres- 
ponds to a constant shear stiffness for the web. The second part above 
the linear part shows a gradual increase in the rate of deformation of the 
web with load (i.e., a gradual decrease in the shear stiffness of the web). 
Since, for all specimens tested, the stress at which the second part of the 
curve started was well below the yield stress for the web material in 
shear, it is reasonable to assume that the change in shear stiffness of the 
web was due to buckling. Since, however, the shear stiffness changes very 
gradually, it is difficult to select consistent buckling loads from plots 
such as figure 4; and in order that the selection be confined to a 
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reasonably short range of load, the effect of a change in the shear 
stiffness of the wet was accentuated in the following manner. 

Method for defining experimental critical stres s.— For each specimen 
tested, the deviation of the load— deformation curve from an extension of 
the linear part of that curve was computed. This quantity was squared and 
plotted against the load on the specimen. The resulting plots of load 
against deviation squared, arranged in the order of increasing stiffness 
ratio 7, are shown in figure 5(a) for one— stiffener specimens and in 
figure 5 (h) for two-stiffener specimens. Each of the curves of the figure 
exhibits a knee which starts as soon as the deformation of the specimen is 
no longer proportional to the load. The buckling load corresponds to some 
point on the knee of the curve of load against deviation squared. 

Comparison of theory and experiment .— A theoretical critical buckling 
load P cr , based on the buckling coefficient from figure 3 and on the 

dimensions listed in table 1, is marked on each test curve of figures 5(a) 
and 5(h) in order to provide a means of direct comparison between test 
data and theory. The subscripts s and c are used to denote whether 
the edges are simply supported or clamped. It will be noted that for l8 
of the 20 specimens the start of the knee of the test curve lies within 
or very close to one edge of the range bracketed by the two extreme 


Maximum s tres ses.— The maximum load sustained by each specimen is 
marked on the curves of figures 5 (a) and 5 (h). For both one— stiffener 
and two— stiffener groups of specimens, the load at failure tended to 
increase slightly with increase in the stiffness ratio 7. In all cases, 
failure ultimately occurred by twisting and collapse of the angle across 
the top edge of the specimen. Since this top angle was of one size for 
all one— stiffener specimens and of another size for all two— stiffener 
specimens, the size and proportions of the longitudinal stiffeners must 
have affected somewhat the maximum load carried by the specimen by 
restraining the top angle from twisting. 


The theoretical shear buckling coefficient for a long flat plate 
reinforced with any number of longitudinal stiffeners, of equal stiffness 
and equally spaced across the plate, can be obtained from a single curve 
for each of the edge conditions — simply supported or clamped. The test 
results were found to be in fair agreement with the theoretical studies. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., December 4 , 19^7 
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APPENDIX 


THEOEETICAL ANALYSIS 


Two solutions applicable to the problem of the shear buckling 
of a flat plate reinforced by longitudinal stiffeners are presented 
in references 1 and 2. The first of these two solutions (reference 1, 
p. 36O) is an approxi m ate solution for the case of simply supported 
edges and is obtained by the energy method. The deflection function 
used was limited to a single half sine wave across the width of the 
plate and did not completely satisfy the conditions of simply 
supported edges along the length of the plate. The solution con- 
sequently yields buckling loads which are too high and unconservative. 
The second of the two solutions (reference 2) presents results for 
the shear buckling of a long orthogonal— anisotropic (often called 
orthotropic) flat plate with either simply supported or clamped 
edges. These results were obtained by solving the differential 
equations of equilibrium of a slightly deflected plate element. It 
is reasonable to expect that this second solution would be applicable 
to the case of a plate reinforced by numerous closely spaced and 
uniformly spaced longitudinal stiffeners. 

It was deemed desirable to obtain a more exact solution than 
given in reference 1 for a long flat plate reinforced by only a 
few longitudinal stiffeners and to obtain some idea of the extent to 
which the solution for an orthotropic plate in reference 2 is 
applicable to a plate with a finite number of longitudinal stiffeners. 
Two energy solutions were therefore obtained. The first solution was 
for a plate with a few or a finite number of stiffeners and the 
second solution, for a plate with a very large number of identical 
closely spaced stiffeners, (in both cases the stiffeners were 
assumed to have some flexural stiffness but zero torsional stiffness.) 

Two different edge conditions were investigated. In the case of 
simply supported edges, the infinite series type of deflection function 
introduced by Kromm (reference 3) was used. This function not only 
provides simple support along the edges but also provides a complete 
set of functions which describe the shape of the deflected surface 
at any section across the plate. Also, in the case of clamped edges, 
a complete set was used. With either function, it is possible to 
approach as closely as desired the exact solution to the problem. 
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SYMBOLS 

h 

width of plate 

°i 

distance from x-axis to ith stiffener 

d 

stiffener spacing 

•N 

•H 

stiffener under consideration 

k 

( (M b2 ) 

shear buckling coefficient \ cr y 

\ « 2 d / 

n, p 

integral number of half waves across plate 

t 

thickness of plate 

w 

deflection in z— direction 

W P 

deflection of plate in z— direction 

w s 

deflection of stiffener in z— direction 

x, y, z 

coordinate axes 

Am.> Bm 

^n* ®n y 
Ap* Bp 

parameters used in deflection function of plate 

D 

/ E t 3 \ 

flexural stiffness per unit width of plate [ p 1 

\ 12 (1 - (i 2 )/ 

E P 

Young’s modulus for plate material 

El 

bending flexibility of stiffener 

Eili 

bending flexibility of the ith stiffener 

F xy 

resultant shear force per unit length acting in middle 
plane of plate 

J 

number of bays across plate 

N 

number of stiffeners 
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T 

work done .by applied shear 

V P 

energy stored in buckled plate 

V s 

energy stored in bent stiffeners 

P 

ratio of half wave length in x— direction to width of 

plate (f) 

7 

ratio of stiffness of single stiffener to stiffness 
of a strip of plate of width d when all ribs are 

identical and equally spaced 

Vi 

ratio of stiffness of ith stiffener to stiffness of 

(tt) 

5 0n 

Kronecker delta (l if n = 0; 0 if n ^ 0) 

® mp 

symbol representing the sum of a trigonometric series 
which takes the value 1 if -- -- is even and if 

— ^ is not even, —1 if m - + ' £ is even and if 
J J" 

m — T ^ is not even^ and 0 for all other cases 

J 

*i> r ± 

parameters used in deflection function of ith stiffener 

X 

half wave length in x— direction 

M- 

Poisson’s ratio for plate material (taken as 0.3) 

V, 6, 0 ± , *i 

Lagrangian multipliers 

F = 2/*^ * 

2 V/ 2 2 


00 



G = m2R m 
m=l 

00 

Hi = — + (— cos 

2 m=l 


b 
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= (l + m2 P 2 ) 2 

K m » = (l + 4m 2 3 2 )' 


Li = __ mt-l) 1 ^ sin — 

m=l 


Wrap = |kP 3 ^_5P — = 0 when m ± p is even 


N 


N, 


mp 


E * 


UUtC-t prtC-t 

sin — - — sin =■ 


i— 1^ 2^ 3^ • • • 


V — / \ rn 2mxcc -* 

Pi = ]T (-1)% sin — ^ 

m=l 


Q = 


_ 


2 P 


= 43 3 k 


■*m 


V = 


(Km*) 2 - (mQ) 2 

Km + 7 

(K m + - (*QV 


R 


00 

-§ + E 




L m 


m=l 


- 


mQ 


(v ) 2 - (W 


V - 


mQ. 


(K m + 7) 2 - (mQ) : 


S = )_ mS m 
m=l 
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Ti = ) m ^ -l)2!IS m COS 

m=l 


2mitcj 

~ 


r — 2mjtc^ 2mnc ^ 

w ij = )_ ** 008 ~V~ sln ~T~ 


m=l 


L iJ 


w 

■L 


2mjtCi # 2mjtc j 


m=l 


R™ sin — - — - sin 

D D 


z ij = 2 + }_ ^ 008 ~ 


2mitC4 2mnc . 


cos 


m=l 


7 

7 = — 

0 J 


V’ = 


2F 


6 * = 


2F 


0. ' = ^ 

2F 


V 

1 2F 


Subscript: 


cr 


critical 
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SIMPLY SUPPORTED EDGES 
Finite Number of Stiffeners 


The deflection function 


v =* sin ^ A m sin + cos 
m=l 



sin 


mjty 

b 


U) 


is assumed to express the buckled shape of the plate. (The coordinate 
system is shown in fig. 6.) The integer m in the deflection function 
represents the number of. half waves across the plate, and the 
parameters A m and B m are associated with the amplitude of the 

mth wave. The width of the plate is b and the half wave length of 
the buckle in the x-direction is X . 

In order to find the energy stored in the buckled plate Vp, 
the energy stored in the bent stiffeners V s , and the work T done 
by the applied shear, the following equations were used: 




d 2 w d 2 w 
dx 2 dy 2 



(2a) 


which is equation (199) of reference 1, and where D is the flexural 
stiffness per unit width of plate; 


V 


s 



dx 


(2b) 


where E^I^ is the flexural rigidity of the ith stiffener and Cj is 

the distance of the ith stiffener from the edge of the plate (only 
bending energy of the stiffeners is considered and the summation is 
extended over all the stiffeners on the plate); and 
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T 



dw dw , 


(2c) 


which is obtained from equation (201) of reference 1, and where 

is the resultant shear force per unit length acting along the longi- 
tudinal edge of the plate . 

It was necessary to consider only the energy and work over one 
half wave length along the length of the plate since the variation 
in deflection is sinusoidal in that direction. 


The deflection function was substituted in equations (2) for 
energy and work, and the indicated integrations were performed with 
the following results : 




A2 N 


(3a) 


m=l 


v 8 - 2 nr 2 

1=1 


00 00 

ZKv, + B m B p) 9ln 


mnc* prtc> 

1 8ln ( 3 b) 


jn=l p=l 


00 00 


T = -2F^ n 


m=l p=l 


m - p c 


(3c) 


where a value for T exists only when m ± p is odd. 

When the buckling load is reached, the structure is in neutral 
equilibrium and is capable of maintaining either the flat or buckled 
form. Mathematically, this neutral equilibrium can be expressed by 
setting the work done by the external load in deforming the plate 
equal to the sum of the energies stored in the buckled plate and in 
the bent stiffeners; that is. 


T « v p + v s 


From this equation, the following critical shear force per unit length 
of the plate is obtained: 
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fv \ * b 3 

V^cr *2 5 . 3 


t 2 5 x 3 


00 00 


mp 


Am ®p 2 2 

. _ m c — -p c 

m=l p=l F 


Z(% 2 *B ^)( 1 + 

m=l \ 


m 




N IT T 00 00 

Z E-Hi's- — m«o 1 pnc. 

2 “5TZ_ Z_ (Vp + b ^p) ain — sln — 


i=l m=l p=l 


(Ua) 


Considerable simplification in the form of equation (Ua) can he 
made by use of the following substitutions: 


P=X b 

%Ii 

7 i = — 


k = 


( F *y) 


cr 


A> 




(l + m 2 3 2 ) 


JL 

iTi 


N mp = 2 ^^ s i n b 


mjtc. pnc. 


sin *- = N 


pm 


Equation (Ua) now reduces to 


00 00 


k = - 


y (*m 2 + X (Vp * B„B p ) R 

m=l p=l 


n m=l 


mp 


83- 


(Ub) 


V, 


mp 


P 2 2 

m=l p=l * “ 
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The wave pattern which causes the shear buckling coefficient k to 
he a minimum is obtained by differentiating equation (kb) with respect 
to each of the parameters A m and B m in turn and setting each of 
the derivatives equal to zero. Two sets of linear algebraic equations 
in Ajjj and B m result. Thus 


*21 (vp + Vp) - 0 (5a) 

p=i 


B *. ("Wp ~ Vp) ’ 0 <5b) 

P=1 


where 


8k(3 3 nrp 
= * m 2 - p 2 



and Mjjp = 0 when m ± p is even. 

Neither equation (5a) nor equation (5b) contains a constant term; 
therefore, in order that the parameters and B m have values 

different from zero (i.e.j the plate takes a form other than flat), 
the determinant in the coefficients of the A m l s and B m *s must be 

equal to zero. The complete determinant is infinite in extent and 
may be represented in the following manner: 


H 

-F=- 


Equation ( 5 &), 

m=l 

A 1 

Ki + N u 

b 2 

M 12 

A 3 

n 13 

B4 

. . 

. . B-l 

, . 0 

a 2 

N i2 

b 3 

0 

A4 

»14 

• • • 

Equation ( 5 b), 

m=2 

-«21 

Kp+ 1$22 

-M23 

N 2 4 . 

»21 

0 

N 23 

0 

• • • 

Equation ( 5 a), 

m =3 

N 31 

m 32 

K 3 +N 33 

M34 . 

. 0 

N 32 

0 

N34 

• • • 

Equation ( 5 b), 

m =4 

-Mui 

n 42 

-M U3 

k 4 +n 44 • 

. Nu 

0 

N43 

0 

• • • 

• 


• 

• 

• 

• 

• 

• 

• 

• 

• 

• 


• 

• 

• 

• 

• 

• 

• 

• 

• 

• 


• 

• 

• • 

• 

• 

• 

• 

• 

• 

Equation ( 5 b), 

m=l 

0 

»12 

0 

n 14 • 

. K1+H11 

-M i2 

N 13 

-M14 

• • • 

Equation ( 5 a), 

m=2 

n 2 i 

0 

N 23 

0 

• **21 

K2+N22 

M 23 

N 2 4 

• • • 

Equation ( 5 b), 

m =3 

0 

N32 

0 

N34 . 

. n 31 

_M 32 

k 3 +n 33 

_m 34 

• • • 

Equation ( 5 a), 

m =4 

N 4 l 

0 


0 

• Mia 

N42 

M43 

k 4 +n 44 

• • • 

• 


• 

• 

# 

• 

• 

• 

• 

• 

• 

• 


• 

• 

• 

• 

• 

• 

• 

• 

• 

• 


• 

• 

• 

• 

• 

• 

• 

• 

• 


= 0 (5c) 


General analysis of determinant Determinant ( 5 c) may be factored as follows: First, multiply 

all the £>dd columns in the right half of the determinant by \1 — 1 and all the even columns in the 
right half of the determinant by —\f— 1, and add each column to the corresponding column in the 
left half of the determinant; and, second, multiply the odd rows in the upper half of the resulting 
determinant by — V^l and the even rows in the upper half by \j — 1, and add each row to the 
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corresponding row in the lower half of the determinant. The lower left quadrant of this altered 
determinant contains only zeros; the upper left and lower right quadrants are equivalent and are 
factors of the original determinant. The factored determinant (infinite in extent) can he 
represented as follows: 


Ki + N U 

m 12 - H 12 

N13 

M lb -\T-* N lk 

m 12 + fZ n 12 

+ ^22 

“ M2 ^ ^ 1 N2 ^ 

N 2h 

n 13 

-Maj -{Zs 23 

k 3 + n 33 

m 3U -\FTn 3U 

*u "it 

»2U 

H 31i + \pT » 3k 

K k + 


It is not readily apparent how the factored determinant might yield a solution in series form; 
but in order to obtain approximations to the true answer for the problem, finite subdeterminants 
can be used . 

A first approximation is obtained by considering only the terms common to the first two rows 
and first two columns of determinant (5d) — the determinant which results from summing m and p 
in equations (5a) and (5b) over 1 and 2. Such a procedure is equivalent to the use of a limited 
deflection function, which is a combination of only one and two half waves across the plate. A 
second approximation is obtained by summing m and p over 1 , 2, and 3, which is equivalent 
to adding to the deflection function used for the first approximation, a term containing three 
half waves across the plate. Similarly, approximations of higher order are obtained by adding to 
the deflection function, terms containing more half waves across the plate. 

ui 
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For example, from determinant ( 5 &) there is obtained for the 
first approximation the following suhdeterminant : 


K 1 + N 11 M 12 - Fl N 12 

M 12 + F 1 N 12 K 2 + R 22 


= 0 


(5e) 


The expansion of determinant (5e) is 


( K 1 + N ll)( K 2 + K 2g) ~ n 12 2 ~ m 12 2 = 0 


( 6 ) 


Equation (6) represents a stability criterion for the shear 
buckling of a long plate with any number, spacing, and size of stiffeners 
along the plate . If the number, spacing, and size of stiffeners are 
known, the proper values of may be substituted in equation (6) 

and the resulting equation solved for k or 7. 

Special analysis of determinant for identical equally spaced 
stiffeners Generally, the stiffeners on a plate are of uniform 
size and are uniformly spaced. Some simplification in the foregoing 
solution may then be Introduced. 

In general. 


If the stiffness ratio 7 ^ is the same for all stiffeners and J is 
the number of bays across the plate, then 


N 



i=l 


J 



i=0 



Dd “P 



J 


mp 
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where 



c J T 4 » El — p . . .. Ill + p , 

o = 1 if — ^ - is even and if — =-^- is not even 
J J 

6 J __ = —1 if — — £ is even and if ™ ~ P is not even 
m P J J 

5^ = 0 in all other cases 

up 


(See reference 1± for derivation of this sum.) 

The determinant (5c) may now be rewritten in terms of 6^" 
rather than N m p. The upper right and lower left quadrants will be 
composed entirely of zeros since das no value if m ± p is 

odd. The upper left and lower right quadrants then remain as 
equivalent factors. The rewritten upper left quadrant (one factor) 
of determinant (5c) appears as follows: 


K i + 

M 12 

r8J 13 

% 

CVl 
1 — 1 
2 

K 2 + 7 5 22 

-“23 

rs j 2 k 

rS J i 3 

-*23 

K 3 + rS J 33 


M U 

7g,J 2l* 

M 3H 

K h + 76 ' 


0 (7) 


S5nce the M^— terms containing the shear buckling coefficient k 

appear in each column, the order in k of the expansion of deter- 
minant (7) is the same as the order of the determinant itself. 


of 


A study of 5 J rr p and determinant (7) shows that the value 

alternates and recurs at intervals and that all the values 


of 5*^ m p are 0 when either m or p is equal to J. Adding and subtracting appropriate rows 
and columns eliminates the 6 J m p-terms (also 7) from all but the first J - 1 columns . As a 

result, the order in 7 of the expression resulting from the expansion of determinant (7) 
is J — 1. It is therefore more expedient to solve this expression for 7 than for k. 


As in the case of the more general determinant (5c), the finite subdeterminants of 
determinant (7) can be used to obtain approximations to the true answer of the problem. Criterions 
resulting from first, second, and third approximations are given as follows for the web with one 
stiffener (J = 2) and with two stiffeners (J = 3): 


First approximation, J = 2 


Second approximation, J = 2 




Third approximation, J = 2 


7 = 


M ^U 2 , M 14 2 , ^ , M 12 2 ( M 12 M ^ * M lU^V 

K 1 K 2 K 3 K 4 


.ffu 


K l K k 


K 2 K 3 


K 1 K 2 


- 1 


1 , 1 (mu - Mu) 2 - Mjg) : 

KjK^ K 1 K 2 K 3 


K 


K- 


(8a) 


(8b) 


( 8 c) 


00 
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«H| CVj 


First approximation, J = 3 


Second approximation 


Third approximation, 



- -jKjKg + K X K U 


7 =i 
' 2 


\/( K l _K 2) 


2 + i+m 12 2 


~ ( K 1 + *2) 


, J = 3 


7 =2 


[ K 1 -*2 + ~YT ) + 4M 12 


K 1 + K 2 " ^ 


J = 3 


[ 23 + 




K^2 + — 


^ K + M 3^ 


- ( m 12 -«iu) ; 


(9a) 


5 

£ 


o 


H 

VJ1 

00 

VO 


+ U 


K 2 m 14 + K 4 M 12 ~ ^ (“23 + M 3t) (^ M 12 M 34 + M lU M 23) 


+ “ 


K 


K7 (”23 + ”3*0 2 + K^3^ 2 + K^3 2 j “ ( M 12 “ M Ht)‘ 


(9D) 


m 
z 

0 12 

1 

tft £ rn 

? o 

7} m <*» 

5 < Z 

n O {,) 

o c 

3§ o 

i/2 - -< 5 
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r * 

Zl *3 
£ 
.-S 
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In equations (8) to (9) the value of 7 is a function of k 
and 0. For any of these equations, y can be computed for various 
vs-lues of 0 at a prescribed value of k. The maximum value of 7 
for the prescribed k can be obtained from a plot of 7 against 0 . 
Curves of k plotted against 7, as shown in figures 3 and 7, are 
obtained by repeating the procedure for other values of k. It 
appears from figure 7 that the third- approximation solution is 
accurate enough and the solution of the determinant of higher orders 
Is not needed. 


If all. the stiffeners are identical and are equally spaced , the 
condition of uniform longitudinal stiffness at any point across the 
width of the plate is approached as the number of stiffeners is 
increased. Consequently, for a very large number of stiffeners, the 
bending energy stored in all the stiffeners may be expressed as 


The work done by the applied shear T and the energy stored in the 
plate Vp are given by the same expressions used in the development 

for a finite number of stiffeners . Following the same procedure as 
that used for the case of a finite number of stiffeners — that is, 
equating energy and work, solving for the critical load, setting the 
derivatives with respect to each of the unknown parameters Aj^ 

and Bjh equal to zero, setting a finite determinant in these 
parameters equal to zero, and solving in this case for the shear 
buckling coefficient k - results In the following equation: 


Equation (10) was obtained by a second approximation, that is, 
by limiting the deflection function to one, two, and three half waves 
across the width of the plate. By a graphical procedure, as previously 
presented, there is found for each value of 7, chosen in equation (10), 
a value of 0 which will make k a minimum. A curve of k plotted 
against 7 obtained from equation (10) is shown in figure 8. For 


Large Number of Stiffeners 




( 10 ) 
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purposes of comparison, there is also plotted in the same figure the 
results of the exact solution for an orthogonal plate from reference 2. 
The second approximation is sufficient to give excellent agreement 
with results obtained from the exact solution. 

Although equation (10) is obtained for the case of a very large 
number of stiffeners, it actually represents the second -approximation 
solution for all cases of three or more stiffeners. This fact can be 
seen from equation ( 7 ) in that the second-approximation solution of 
equation ( 7 ) for all cases of three or more stiffeners gives the same 
results as equation (10). 


CLAMPED EDGES 

Finite Number of Stiffeners 


The problem is more involved in the case of clamped edges. A 
different deflection function is assumed and the Lagrangian multiplier 
method (references 5 and 6) is used. With the new coordinate system 
shown in figure 9> the deflection function for the plate is assumed 
to be expressed by 


w p = sin ^ sin + C os B n cos (11) 

m=l n=0 


This expression is a complete set of functions symmetric with respect 
to the origin. Since the plate is infinitely long, the expression of 
deflection by a complete set of anti symmetrical functions will give 
the same results. In order to ensure zero deflection and zero slope 
at the edges, the expression Wp (equation (11)) is subject to the 

following restraining conditions: 


00 

X 

n=0 


(“if ^ = 0 


(12a) 


4 “if % 


m=l 


0 


(12b) 
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A set of new deflection functions is used for the longitudinal 
stiffeners. For the ith stiffener. 


( w .) 1 * *i 


sin 2S. + r 


C08 


nx 

X 


(i = 1,2, ...N) 


where N is number of stiffeners. In order that the deflection 
the plate directly under the stiffener and the deflection of the 
stiffener itself be the same, the expressions for w p and ^w a ^ 

are subject to the following restraining conditions: 


Z 2mjtc-» 

V sin — = 0 (i = 1,2,... N) 

m=l 


r 


I 



cos 


2nitC4 

i- = 0 

b 


(i = 1,2,... N) 


where c^ is the distance from the x-axis to the ith stiffener. 

When the expressions for w p and w 8 are substituted in 
equations (2), the energy expressions become 



00 


00 


X VN , 2 * Z VV 2 



v a . 2F 



i=l 



T 


F xy * \/_ 
m=l 


(13) 

of 


( l^a) 


(l^b) 


( 15 a) 


(15b) 


( 15c ) 
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where 


F 


= & / n\ 4 X h 

2 W 2 2 


V = (1 + ^p 2 ) 


2 


3 


x 

b 


?i 


E l T i 

Dd 


On - 1 

if 

n = 0 

O 

II 

£ 

if 

n 0 


The total energy V = Vp + V 8 — T is then to be minimized. During 

the minimizing process the restraining conditions, equations (12) 
and (lh), can be satisfied by the Lagrangian multiplier method. (See 
reference 5») 

The following notation is used for the Lagrangian multipliers: 
v and € correspond to equations (12a) and (12b), respectively; 
and correspond to equations (lUa) and (14b), respectively. 

Then the function to be minimized is 


w W 

f = Vp + v 3 - T - (-if B n - 6^2 (-!f “Am 


n=0 


m=l 



(16) 


If the function f is minimized with respect to each of the 
parameters A^ B n , Aj, and Tj in turn, the following expressions 

can be obtained: 


2k 
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df 


N 


2 . m V~ , 2mnc> 

= 2FK m , A m - F xy n mB m - era(-l) + y 0j_ sin — - — 

= 0 ( m = 1, 2, ... 00 ) 


(IT) 


af_ 

3b 


n 


- Z^n'C 1 * & 0n) B n - »xy« 2 «*n 


N 

> 

1=1 


Z 2nnc, 

*i 008 — b = ° ( n = 0,1,2, ... 00 ) (18) 


ij- - K) - «i 

= 0 (i = 1,2, ...N) ( 19 ) 

% - ^iK) - *i 

= 0 (i = 1,2, ...N) (20) 


Equation ( 18) can "be separated into the following parts : 

N 


UFB 0 - v + 


]>_ *1 » 0 (c 


orresponding to n = 0) (21a) 


i=l 


2FVB m - Fxyn 2 ^ - V(_i) m + cos = 0 (21b) 


i=l 


(m = 1,2, . . .“ ) 


Solving equations (IT) and (21) together and using the notations V* = — , 

€ * = and 'l' , ' = result in the following expressions 

2r 2F 2F 

for the A*s and B's: 


WACA TW Wo. 1589 


25 


% = R m 


+ S 


l) m - sin 


m 


1=1 
N 


V(-l) m -^T 


• cos — -i 


1=1 


(m = 1 , 2 ,. ) (22) 


B m = 
in m 


+ s 


N 


2mnc, 


C08 


1=1 


m 


N _ 

, m(-l) m -^' 0j» sin — - 1 
1^1 


2uutc J 


( in = 1,2, ...<*>) (23a) 


N 

B ° *f v ‘ -5X* 1 ' 


1=1 


(23b) 


where 


R m = 
m 


V 


(v ) 2 - 


Sm = 


mQ 


(Km*) 2 - (*Q)‘ 


q = = 1+3 3^ 

2F 


If equations (19) and (20) are substituted In equations (ll), the 
following expressions result: 


1 

27, 


2mnc* 

-h' -/ Am sin = 0 (1-1,2,...N) 

1 m=l 


(21a) 


00 

B n 


2n«c J 


cos 


= 0 (1 = 1,2, ...N) (2Ub) 


n=0 
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If Ajjj and B n as given By equations (22) and ( 23 ) are put into 

equations ( 12 ) and (24), the following set of simultaneous equations 
is obtained with v*, 0^*, and \|rj* as the variables: 

N N 

v «r + € »s + ]T h' (~ p i) + Yi ** (“*0 = 0 (25a) 

1=1 fe 

N N 

V *S + e»G +^_ 0 1 *(-l 1 ) + ]T V(- t i) - 0 ( 25 b) 

1*1 i=l 

H N 

v’(-Fj) + « * (-*- j) ♦H*l , ( i «) + **J* 577*^1 ♦i’( w u) ' 0 <25=) 

i=l J i=l 

(J = 1,2,... N) 

N N 

+ ,'(-Tj) -^X^) + ♦j'fir - ') * 0 ^ 

i=l i=l v 'J'' 

( J = 1,2,... N) 


where 


B 


1 

2 


00 



w 

■Z" £ 


m 




00 

Hj - f + ]T 


2mjtc j 


m cos — 


m=l 


L J " 


„ 2mnc , 

m(-lfB m sin 
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2nufc 2nufc 

[ u =2_ R * 8in sin V = x 

m=l 


Ji 


2mnc 2mnc 

1 --- J. = z, 


J = I + y R m 008 — T“ 008 ~b~ = "ji 


m=l 


00 

-z 


n£ 


m 


m=l 


Z rri 2 mrtC 4 

(-1) S„ sin 

m=l 


w 

T j = 2Z 008 


2mjic , 


m=l 


00 2mitc, 2mnc , 

W iJ = 2 _ S m C ° 8 T“ 8in "b 
m=l 


Thera are 2N + 2 equations in the preceding system of simultaneous 
equations and there are 2N + 2 unknowns, namely, €*, , 

^ 2 ', **• ^n'* Xn order to ensure a 

nonvanishing solution for these Lagrangian multipliers, the following 
determinant must he zero: 
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ro 

00 

Equation 

( 25 a) 
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s 
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- P 2 
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( 25 b) 
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- L 1 

~ T 1 
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. . . -Lj, 
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Equation 

(25c), 

J=1 
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W 11 

X 21 
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... x N1 

W N1 


Equation 

(254), 

J=1 
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- T 1 

W 11 

Z*l *1 + ■ ■ 
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W 12 

Z 21 

... . 

Z N1 


Equation 

(25c), 
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" P 2 

-l 2 

X 12 
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w 22 

• 

... x N2 

W N2 

s 0 ( 26 ) 

Equation 

(254), 
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# 
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• 
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• 

• 
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• • 
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W NN 
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(254), 
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W N2 

Z 2H 

. . . w jjn 
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Determinant (26) Is of the 2N + 2 order; it can be greatly simplified by using the numbering 
system shown in figure 9 for the longitudinal stiffeners, which are assumed to have equal stiffness y 0 
and are equally spaced across the width of the plate. In this system, the first stiffener, if any, 

1 8 the central stiffener; the numbers 2, 4, . . . denote the stiffeners on the positive y-eide of 
the plate, and the numbers 3 * 5 » • • • indicate the stiffeners on the negative y-eide, numbering 
out from the center in each case. With this numbering system and where p = 1, 2, 3> • • •> the 
following relations follow immediately: 
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c 2p 

= ~ °2p+l 

Z 11 

= R 

H 2p 

= H 2p+1 

ro 

= “ P 2p+1 

T ?P 

= T 2p+1 


= _L 2p+l 

2p,i 

= W 2p+l,i 

X 2p,i 

= _X 2p+l,i 

'2p,i 

= Z 2p+l,i 

W i,2p 

= -*1,2?+! 


With these relations, the determinant can "be simplified by means of operations of the 
following type: First, columns under 0 2 p* are &d-<t©d- to columns under 0 2 p+i* > then rows 

corresponding to equation (25c) with J = 2p + 1 are subtracted from rows under equation (25c) 
with J = 2p. In this way, all the columns under ^op+q* are eliminated in determinant (26). 

Similar operations can be used to eliminate all columns under ^pp + q' • After dividing some 

rows for convenience by 2, the resulting determinant becomes 
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R 
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~ p 2 

-h 2 
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O 
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~ P 4 

- l h 
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*24 

'O 
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-TU 

z lU 

*42 

Z 2U 

*44 

z '* u + >7 


= 0 (27) 


Determinant (27) Is now of the N + 2 order. For the particular cases 
of one and two longitudinal stiffeners, the corresponding determinants 
are as follows: 

For one stiffener 


For two stiffeners 
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-T x R + 


1 

?7o 


= 0 


( 28 ) 


R 

S 

CM 

(X. 

1 

CM 

W 

1 

S 

G 

CM 

1 

1 

^3 

ro 

“ P 2 

-Lr 

X P ? + r^~ 
22 ^0 

*22 

-h 2 

- T 2 

CM 

CM 

z 22 + ir 


'+7o 


= 0 


(29) 


MCA TN No. 1589 


31 


In determinants (27) to (29), each element is an infinite series 
consisting of functions of k and 3. For each of the determinants, 
7 0 can he computed for various values of 3 at a prescribed value 
of k. A maximum value of y Q for this prescribed k can be 
obtained by plotting y Q against 3 . A number of these maximum 
values of y Q corresponding to various values of this prescribed k 
can be obtained. Instead of plotting k against y Qt k was 

plotted against 7, as was done for the case of simply supported 
edges, where 7 = J7 Q and J equals the number of bays into which 

the plate is divided by the stiffeners. Thus, J = 2 for one 
stiffener and J = 3 for two stiffeners. The advantage of such a 
plot is that the various curves for different numbers of stiffeners 
become almost coincident. (See fig. 3.) 


Large Number of Stiffeners 

In the case of a very large number of stiffeners of equal 
stiffness and equally spaced across the width of the plate, the 
bending energy stored in the stiffeners may be expressed. Just as 
in the case of simply supported edges, as follows I 


V 


s 


P b/2 

1 El / 

2 d / 

J- b/2 




dx dy 


If, now, the same deflection function as given by equation (11) for 
the plate is used for the stiffeners, the following expression is 
obtained for the bending energy stored in the stiffeners: 


V 


s 



00 


+ X B n 2 ( X + 
n=0 



( 30 ) 


There is no change of expressions in Vp and T. Of course, since a 

separate deflection function for the stiffeners is no longer used, 
there is no need for the restraining conditions of equations (14). 

The function to be minimized is 


f = V„ + V a 


_ T - v ^ (-if B n - € m (-DX 


n =0 


m=l 


( 31 ) 
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If the sane steps are followed as before, the following determinantal 
equation is found to be the condition for the existence of nonvanishing 
solutions for the Lagrangian multipliers: 


y~ 

m=l 
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mS » 


m=l 


'"6m 1 


m=l 


= 0 


m 2(1 + 7) 


V 


m=l 


(32) 


where 
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V + r 

(V + y ) 2 - M ) 2 


s 
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mQ 

(v + r ) 2 - (mQ) 2 


In determinant (32), each element contains an infinite series of 
functions of k, 7, and 3. For prescribed values of both 7 and 3, 
several values are assigned for k and the corresponding values of the 
determinant can be determined. If the values of the determinant are 
plotted against k, one value of k can be found which makes the value 
of the determinant vanish. This particular value of k is called ko . 
Now, if the value of 3 is changed (while the value of 7 remains 
the same), the corresponding value of kg is also changed. There 

exists a certain 3 which makes the value of ko a minimum. This 
minimum value of k Q is the critical shear buckling coefficient k 

corresponding to the prescribed value of 7. In a similar manner, 
other critical shear buckling coefficients can be determined for other 
prescribed values of 7 . Finally, a curve can be obtained with k 
plotted against 7. This curve is presented in figures 3 and 8. In 
figure 8, the exact solution for an orthotropic plate from reference 2 
was also plotted for comparison. 

From the results obtained in the case of simply supported edges, 
it is believed that the curve obtained for a large number of identical 
equally spaced stiffeners represents the solution for all cases of 
three or more stiffeners. 
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DISCUSSION OF THEORETICAL RESULTS 


In the foregoing analysis it is seen that for both simply supported 
and clamped edges a determinant is obtained for the determination of the 
shear buckling coefficient k. In the case of simply supported edges, 
the usual energy method was used and the determinant has been solved by 
first, second, and third approximations. They converge very fast, 
however, as may be Judged from the results of the solution for one 
stiffener. (See fig. T») Also, figure 8 shows that the second-approxi- 
mation solution for a large number of stiffeners is sufficient to give 
excellent agreement with results obtained from reference 2. In the 
case of clamped edges the Lagrangian multiplier method was used and 
the solution of the determinant is exact. If enough terms In each 
infinite series are taken, the solution can be made to any desired 
accuracy. 

Curves of the stiffness factor 7 plotted against the shear 
buckling coefficient k for one centrally located stiffener, two 
Identical equally spaced stiffeners, and a large number of identical 
equally spaced stiffeners are shown in figure 3 for both simply 
supported edges and clamped edges . The curve for a large number of 
identical equally spaced stiffeners represents the solution for all 
cases of three or more stiffeners. Since the curves of figure 3 do 
not depart from one another by more than about 2 or 3 percent over 
the range shown, practicability would dictate the U3e of the curve 
for a large number of identical equally spaced stiffeners to predict 
the buckling load for a plate with any number of stiffeners, provided 
the buckling coefficient so obtained is not higher than could be 
obtained by replacing the stiffeners by simple supports. On the basis 
of the criterion Just stated, it is apparent from figure 3 that in 
the case of a simply supported plate reinforced with a single stiffener 
no further increase in the shear buckling coefficient k can be 
obtained by increasing 7 beyond about 2000. Similarly, it Is seen 
that for the case of a simply supported plate reinforced by two 
identical equally spaced stiffeners no further increase in the shear 
buckling coefficient k can be obtained by increasing 7 beyond 
about 14-5,000. 
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TABLE 1 

DIMEWS IONS OF SPECIMENS 


Specimen 

Plate 

width, 

b 

( in.) 

Plate 

thickness, 

t 

(in.) 

Plate 

length, 

(in.) 

Stiff- 

ener 

height 

(in.) 

Stiff- 
ener 
width 
( in.) 

Stiff- 

ener 

thick- 

ness 

(in.) 

7 - S 
7 Dd 



One— 3tiffener 

specimens 



1 

6.00 

0.0340 

48.00 

— 

_____ 

_ _ 

0 

2 

6 .04 

.0334 

47.98 

0.185 

0.31 

0.0281 

16.5 

3 

6.02 

.0337 

48.05 

.300 

• 315 

.0376 

78.4 

4 

6.00 

.0330 

48.12 

.305 

.290 

.0331 

153 

5 

6.01 

.0325 

48.09 

.420 

.460 

.0391 

234 

6 

5.99 

.0334 

48.00 

.510 

.385 

.0494 

477 

7 

6.00 

.0326 

48.07 

• 550 

.500 

.0391 

503 

8 

6.00 

.0327 

48. 06 

.550 

.650 

.0492 

633 


5.96 

.0313 

48.00 

.622 

.642 

.0395 

824 

10 

5.98 

.0316 

48.02 

.581 

.653 

.0626 

1050 

11 

5. 96 

.0315 

48.04 

.623 

.718 

.0621 

1290 

12 

6.00 

.0315 

48.05 

.668 

.700 

.O63O 

1570 

13 

6.00 

.0312 

48.03 

.665 

.715 

.0622 

1590 

14 

6.06 

.0330 

48.02 

.985 

.795 

.0618 

4180 



Two— stiffener 

specimens 



15 

7.88 

0 .0326 

63.O 

— 

— 


0 

16 

7.88 

.0326 

63.O 

0.248 

O.505 

0.0327 

53-4 

17 

7.88 

.0321 

63.O 

.348 

.649 

.0334 

147 

18 

7.88 

.0326 

63.O 

• 550 

•721 

.0328 

484 

19 

7.88 

.0325 

63.O 

.685 

•794 

.0510 

1390 

ao 

7.88 

.0325 

63.O 

.944 

.783 

.0634 

4540 
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Section A- A 

(a) Specimen with one 
stiffener. 


Section C-C 

(b) Specimen with two 
stiffeners. 


Figure 1.- Test specimens. (All dimensions are in inches.) 
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Figure 2.- Photograph of test setup. 
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Shear deformation, in. 

Figure 4.- Typical load -deformation curve for shear webs 
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o Jest results 

(&) (Per) Computed shear buckling load Per damped 
c ’ 5 and simply supported edge condition , respectively 



(a) Webs with one stiffener. 

Figure 5. - Comparison of test results and theory on the shear buckling 
load of long flat plates with longitudinal stiffeners, identical and 
equally spaced. 


k2 
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(b) Webs of two stiffeners. 
Figure 5.- Concluded. 
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1+3 



are indicated) 



Figure 6.- Coordinate system used in theory for simply supported plates. 



Figure 7.- Variation .of shear buckling coefficient with rib stiffness for simply supported 

plate with one longitudinal stiffener on center line. First, second, and third approximations. 
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Figure 8.- Comparison of results of energy solution of present paper with results 

of exact solution from reference 2 . 


4=r 

VJ1 


NACA TN No. 1589 


46 


NACA TN No. 1589 



Figure 9.- Coordinate system used in theory for clamped plates. 


